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ABSTRACT 

Screened potential, modified by non standard electron cloud distributions responsible for the shield- 
ing effect on fusion of reacting nuclei in astrophysical plasmas, is derived. The case of clouds with 
depleted tails in space coordinates is discussed. The modified screened potential is obtained both from 
statistical mechanics arguments based on fluctuations of the inverse of the Debye-Hiickel radius and 
from the solution of a Bernoulli equation used in generalized statistical mechanics. Plots and tables 
useful in evaluating penetration probability at any energy are provided. 
Subject headings: Plasmas — nuclear reactions — atomic processes 



1. INTRODUCTION 

Nuclear fusion reaction cross sections and rates 
are sensitive to the screening effect of the electron 
cloud around reacting nuclei, an effect that has been 
widely investigated both theoretically and expcrimcn- 
tally since the early works of Salpeter (iSalpe tcr 1953; 
iSalpeter fc Van Hornl[l969h . 

Different situations arise when fusion reactions take 
place: 1) in laboratory experiments, where a metal or 
gaseous target of a given element is bombarded by an 
ionic or charged particle beam, electrons are for the most 
part bound in atomic orbits and few of them can be con- 
sidered free; 2) in stellar cores and other space and as- 
trophysical plasmas where ions and nuclei are embed- 
ded in an electronic environment made by mainly free 
electrons; 3) in a deuterated metal and other solid-state 
matrices where an impinging deuteron beam reacts with 
implanted deuterons. 

— In laboratory experiments, penetration through a 
screened Coulomb potential at center of mass energy E 
is shown to be equivalent to that of bare nuclei at en- 
ergy E + Uc where Uc — Zi Z2e'^ /R^ and i?a is the 
atomic radius or the radius of the innermost electrons. 
Zi e and Z2 e are the charges of the two reacting nu- 
clei and Ue is usually taken as constant in the evalu- 
ation of cross sections and rates at any energy. Very 
often, at very low energy, fusion cross sections mea- 
sured in laboratory experiments are higher than the value 
calculated by means of the usual Debye-Hiickel (DH) 
screening factor. Stopping power of the incoming beam 
and temperature reached after energy deposition are im- 
portant quantities for the correct measurement of the 
cross sections. Modifications of electron distribution 
can be induced. The screening effect must be evalu- 
ated to obta in the correct a strophysical factor at very 
low energies (jAssenbaum et al . 1987; Carraro et al.lll988l : 



Bracci et al.] 119901: iShoppa et all 119931 : iStrieder et all 



2OOII) . 



— In astrophysical plasmas, which that can be con- 
sidered ideal, the free electrons move around the react- 
ing nuclei and occupy a sphere of DH radius R^^^ = 
yJkT / [At: e^ nZp) which is taken on the order of i?a, 

with n the particle density and Zp = J^ii^f + ^i) ^i/^it 
where the sum is over all positive ions and Xi is the 
mass fraction of nuclei of type i. Only with a de- 
creasing radius does the screening effect become impor- 
tant. A screening factor of the rate can be derived 
when the energy of the Gamow peak E^ > Uc- It 
is given by the Debye factor / = exp(Uc/kT) with, 
this time = Z^ Z^e^/R^^ ([RoTfs & Rodnev 200i 
Ichimarul [1991 ICastellani et all Il997l : lOpher fc Opherl 



200Cl| ). Recently, it has been clarified through numer 



Electronic address: 
Electronic address: 



piero.quarati@polito.it| 



aSomo!scSfone<aipolito.it| 



ous experimental observations that the velocity distri- 
bution function of electrons (and possibly also of ions) 
in stellar atmospheres, in space and astrophysical plas- 
mas, may deviate from a Maxwell-Boltzmann distribu- 
tion in the high-energy tai l if non-local thermodynamic 
effects are non-neglig i ble dO xcnius' ICoIIi^ fl98l 

lPevraud-Cuencal[l992l : IChevaUier 200 61 . In stellar atmo- 
spheres, atomic processes such as radiative and dielec- 
tronic recombination show rates depending on deviations 
from Maxwellian distributions of electrons (Macr o et al.l 
I2OO6I ). In stellar cores, signals of possible deviations of 
ion distributions are evident and, although small, should 
be considered because they are capab le of sensibly in- 
fluenc ing nuclear fus ion reaction rates (|Ferro &: Quaratil 
2005; Lissia fc Ouar ati 2005). 

— In deuterated metals or solid-state matter, the 
strong screening effect has still to be clearly under- 
stood and discussed, although a few interesting descrip- 
tions have recently b een brought forward to reproduce 
experi m ental r esults (iRaiola et al.l 120041 : ICoraddu et al.l 
|2004a"b', '200i iKim fc Zubarevll2006[ ). The approach we 
are discussing here is very useful in understanding the 
fusion rates in this matter which could simulate some 
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high-density astrophysical plasmas. However, this apph- 
cation deserves a separate detailed paper and we will not 
discuss this case here. 

An important issue in the shielding of electrostatic 
potential in plasmas concerns the investigation of non- 
linear charge screening effects that can induce modifica- 
tions in the DH pote ntial, usually derived by linearizing 
the P oisson equation (jCraven il997L [G ruzmov fc Bahcalj 
|1998() . One of the first studies which, on a micro- 
scopic basis, demonst rated dev iations from DH factor 
can be found in (Johnson et a l.in_992). [See also (.Shaviv, 
I2OOI iShaviv et al l l20nnl . 12001. 2003; iChitanvisI [2007f )]7 
The different approaches elaborated up to now are 
based on the assumption that electrons are distributed 
in space according to a Boltzm ann fac tor. Few au- 
thors (Bryant 1996 : P&euman et al. 2004; lLeubneill2004l : 
iKim fc Juntri2004l : iRubab fc Murtazal I2006D have as- 



sumed a stationary A;-Lorentz distribution to describe 
significant deviations from standard distributions, con- 
structed from experimental distributions, due to the 
prese nce of enhanced high-ene rgy tails. In a recent refer- 
ence (jRubab fc Murtazall2006[ ). for instance, we can find 
quotations from many works where such distributions are 
reported. By analyzing energy profiles the authors derive 
an effective length smaller than the standard DH radius, 
depending on the fc-parameter. This result is obtained 
by solving a Poisson equation in the linear weak approx- 
imation. The consequence must be an increased barrier 
penetration factor. 

Our approach differs. We assume that the electron 
cloud is spatially distributed following a generalized 
steady state distribution of the q-type that reduces to 
an exponential distribution when the q parameter (also 
known as the Ts allis entropic parameter) approache s 
the g ^ 1 limit dTsallid [Hsl iTsaUis & Borges '2003^. 
We refer the reader to Refs . (Lcubncr fc Voros 2005; 
lLeubneill2005l : iBurlaga et al. ''2006) for a detailed descrip- 
tion of Tsallis generalized statistics and some of its ap- 
plications to astrophysical problems. We then calculate 
the modified screening potential by considering two dif- 
ferent approaches. One concerns the use of a general- 
ized Poisson equ a tion o r Bernoulli equation as used by 
iTsallis fc BorgesI (l2003l); th e other is based on super- 
statistics fe eck 2001L 120041 : iWilk fc Wlodarczvkl [200I 
IWilkl [20061 ) considering fiuctuations of an intensive pa- 
rameter (the inverse of the DH radius). This implies fiuc- 
tuation in temperature and density of the components 
of the plasma. Once we obtain the modified potential, 
we calculate penetration probability through that poten- 
tial. In this paper we limit ourselves to the value range 
< q < 1 (distribution in spatial coordinates with de- 
pleted tail cut at about {l/R^^^)), leaving the case q > 1 
that describes a distribution in space coordinates with 
enhanced long-distance. 

In the standard DH shielding approach, after lineariza- 
tion of the Poisson equation, the electrostatic screened 
potential behaves like T/^^jj ~ exp(— r/i?^^) where 
r is the coordinate with which to evaluate the DH po- 
tential. The assumptions made to deriv e the above re- 
lation are, among others things, that (jCravensI Il997l : 
lBellanl[2006l ): plasmas are colhsionless; the induced per- 
turbation is slow and depends slowly on time (slowness); 
only electrostatic fields are present while induced fields 
are negligible; temperature is spatially uniform and the 



plasma remains in equilibrium during a perturbation; a 
temperature can always be defined; the number of par- 
ticles inside the DH sphere is large and therefore fluc- 
tuations are small; although ions and electrons have a 
random thermal motion, perturbations induced around 
the equilibrium, responsible for small spatial variations 
of electrostatic potential, can be neglected. 

Of course, if the real situation differs from the one im- 
posed by one or more of the assumptions reported above, 
the use of DH potential may induce errors in the evalua- 
tion of the penetration factor and nuclear reaction rates. 
Deviations from the conditions imposed by the assump- 
tions are taken into account in this work by choosing the 
inverse of DH length 1/i?^^ a fluctuating parameter. 
The electrostatic quantity r V{r) is asymptotically given, 
in this case, by a power law instead of an exponential 
law because it must satisfy a differential equation, the 
Bernoulli equation (or a special case of it), that has power 
law functions as solutions. 

We can derive the modified DH potential Vq{r) and, a 
posteriori, the charge distribution pq as an asymptoti- 
cally power law function. Also in the standard DH ap- 
proach also two equations are needed, one coming from 
electromagnetism (Poisson equation), the other from sta- 
tistical mechanics (Boltzmann factor). 
Penetration factor T{E) can be calculated by means of 
the WKB approach using the modified DH potential. We 
obtain results that differ from the ones calculated with 
standard DH potential and that will be useful in the 
interpretation of experimental results on atomic and nu- 
clear rates in several astrophysical systems and processes. 
We also derive equivalent energy Uq that we give by an 
interpolating analytical expression, a plot and a tabula- 
tion. Energy Uq results a function of the variable D/E, 
where D is defined by D — Zi Z2 (l/R-^^). It is easy 
to observe that Uq depends, for a wide range of values of 
D/E, on l/^/kT and V^. 

In Section 2, we explain how to derive the modified po- 
tential on the basis of super-statistics arguments with the 
inverse DH length as a fiuctuating parameter; we then in- 
troduce a nonlinear differential equation, to be associated 
with the Poisson equation, whose solution coincides with 
the potential derived directly from the super-statistics 
approach. 

In Section 3, we derive the penetration factor that can 
be used for evaluation of nuclear fusion rates in astro- 
physical plasmas and indicate the range of validity of 
approximations adopted. In Section 4, we discuss some 
representative examples whilst, in Section 5, we report 
our conclusions. 

2. MODIFIED DEBYE-HUCKEL POTENTIAL 

Combining the Gauss law and the relation that links 
the electrostatic field to the electric potential V{r) of a 
point test unitary charge at the origin in a vacuum, from 
the Poisson equation one obtains the pure Coulomb po- 
tential. After a sufficiently long time, electrons and ions 
rearrange themselves as a response to the forces on them. 
Ion density eventually remains uniform while electron 
density near the test charge increases. At the new ther- 
mal equilibrium, the distribution of electrons in an elec- 
trostatic field is assumed to be given by the well-known 
Boltzmann factor. Assuming the Boltzmann factor for 
all the particles, after linearization and using neutrality 



Modified Debye-Hiickel Shielding 



3 



condition, the Poisson equation can be written as 



1 



i?2 

DH 



Vir) , 



(1) 



with solution given by the DH potential and the charge 
density p^^ expressed as 

; 

When one or more of the linearity constraints are vio- 
lated or relaxed, a different description of the screening 
is required. 

If we assume that nonlinear effects produce fluctua- 
tions on the inverse DH radius, by following the approach 
usually develope d by super-statisti cs for inverse temper- 
ature P = l/kT (lBeckll200lL[200l . we can describe the 
plasma around the test charge as made of cells where 
is approximately constant and the system can be 
described by ordinary statistical mechanics, in this case 
by the exponential (Boltzmann) factor exjp(—r/Rj^^). In 
the long term run the system is described by a spatial 
average over the mean of fluctuating quantity 
Fluctuation of the inverse DH radius also means fluctu- 
ation of the plasma parameter given by 

3/2 



1 



7 ^ 



4 TT 
kT 



(3) 



With few changes, we follow the approach by Wilk and 
Wlodarczyk for the case o f distributions with deplete d 
tails {q < 1) dWiig [20061 : iWilk fc Wlodarczvkl [200l . 
Here we focus our attention on g < 1 distribution be- 
cause this shows a depleted tail with a cut-off, that is, 
the spatial distribution we assume for the electrons. 
We assume that a certain variable r of the system is lim- 
ited between and [(1 — q) Ao]~^ where Ao is a constant 
parameter. 

We define the function 

oo 

^q<i{r, Ao) J fq<^ ^5 -^o) cxp (-Ar) d\ , (4) 



where Cq is a normalization factor and fq<i{r, A; Aq) is 
the probability density to observe a certain value A of the 
system which is spread around the value Aq. The expres- 
sion we choose for /g<i(r, A; Aq) is a gamma distribution 

/,<i(r. A; Ao) = ^"^''^ ^ At^-^ exp (-A A,(r; Aq)) 



where 



^q{r; Ao) = ^— ^ Aq ^ - r , 



(5) 
(6) 



and q is the entropic Tsallis parameter. 

Inserting the function /g<i(r. A; Aq) into Tq^i{r; Aq) we 

obtain the normalized power law distribution 

.F,<i(r; Ao)-Ao (^l-^Aory^ . (7) 

Average value and variance of A depend on variable r 
according to 

1 — 2-g 



A = 



{l-q)Aq{r- Ao) ' 



A2 = 



[{l-q)Aq{r- \,)\ 



(8) 



where x{r) — J x{r, X) fq<i{r, A; Ao)(iA is evaluated by 

means of distribution fq<i{r, A; Ao). 

However, the relative variance depends only on q 



(jj 



(A)2 



(9) 



We remark that quantity Ao coincides with the spatial 
average of A, that is 

(A) = Ao , (10) 



where (x) = J x{r) Tq<^i{r; Aq) dr is evaluated by means 
of distribution .Fq<i(r; Ao) which is the weighted average 
of the exponential (or Boltzmann-like) factor exp(— Ar) 
with weight equal to /g<i(r, A; Ao) and coincides with 
the Laplace transform of fq<i{r, A; Ao). 
Some special limiting cases are 



and 



fq^iir, A; Ao) = S{X - Ao) . 
Tq^i{r: Ao) = Ao exp(-Aor) 



(11) 
(12) 



the ordinary Boltzmann-like factor, with (A^) = (A)^ and 

fq^oir, A; Ao) = Ao{r; Ao) exp - XAo{r; Ao)) , (13) 

where Ao^r; Ao) — 2/Ao — r and 

^,^o('^; Ao) = Ao(l-Aor) , (14) 



a linear function, with (A^) — 2 (A)^. 
Let us now identify functional J^g<i(r; Ao) with quan- 
tity rVq{r) and substitute A with l/iJ^j^, so that Ao — 
(l/i?Qjj) coincides with the spatial average of the DH 
radius fluctuation. By posing 



C, = (2-g)(l/i?„J~i 



(15) 



a characteristic length of the system under inspection 
which reduces to R^j^ in the g — > 1 limit, we obtain the 
following expression 



Vqir) 



1 



l-(l-g) 



(16) 



which, for g ^ 1 reduces to the standard DH potential. 
We also have charge distribution 



Pqir) 



1 



r \ 



(17) 



which for 9 — s- 1 reduces to p^^ , the charge distribution 
of the DH approximation. 

Therefore, by considering l/R-^^ subject to fluctuations 
described by a gamma distribution, quantity rV{r) re- 
lated to the potential energy barrier is modified from 
the exponential DH expression to a power-like law, typ- 
ical of generalized q < 1 distribution. When we can 
establish the functional relation between density n and 
temperature k T, as, for instance in a solar-like star, 
whe re quantity n/(fc r)'^ is constant along the star pro- 
file (iRicci et al.l[l99l . by means of the relative variance 
we can establish a link between inverse DH radius and 
temperature fluctuations and parameter q through 



A(l/i?„J _ Ay^VfcT _ A(fcT) 



1/R, 



^/^{JkT 



kT 



= {l-qf'^ . (18) 
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Diffusion of matter between layers with different temper- 
atures induces local temperature fluctuations and den- 
sity perturbations; fluctuations around an equilibrium 
or a steady-state matter profile induce fluctuations of 
quantity in the DH sphere, particularly in the 

regions where the number of particles inside the DH 
sphere is small. Density and temperature fluctuations do 
not alter the macroscopic plasma parameters and must 
agree with the requirements of the constraints imposed 
by the ma croscopic observations. In the solar interior, 
iGruzinovfc Ba hcall [1998.) have evaluated that the DH 
radius is about 2 • 10~^cm, therefore containing a small 
number of particles with a corresponding non-negligible 
particle fluctuation. 

Let us now justify expressions (jl6p and (|17p by means 
of another approach that considers a generalized version 
of the Poisson equation. In fact, in the standard case, 
the DH potential can be obtained from the solution of 
the second order differential equation of the type 



dy 

T ="'y 

ar 



y 



(19) 



with y = rV{r) — exp(ar). 

We replace the above linear equation with the following 
one 



or 



d'y 



dr 



with 



y = exp (agr)= 1 - {1 - q) ag 



(20) 



(21) 



where g is a real number parameter and coincides with 
the Tsallis parameter. For q — s- 1 we obtain exp]^(ar) = 
exp(a r). 

To be explicit, for our case we generalize Eq. ([T]) into 



d^ 

dr'^ 



2q-l 



(22) 



which is a generalized Poisson equation whose solution 
coincides with Eq. (|16p . 

Wc report, for completeness, that by considering 
the Bernoulli equation introduced by iTsallis fc Borgea 

IIOM) 

dy 

= aiy + Ugy"^ , 



dr 

or 

d'^y 

— afy + {1 + q) ai aqy'' + qugy 
the solution is 



2 ,,2 q-l 



y 



ai \ J 



(23) 



(24) 



(25) 



By posing y — rV{r) we obtain 

rV{r) = l + {l-q)ag'i 
when ai = and r V{r) = exp(ai r) when — 0. 



3. PENETRATION PROBABILITY 

We calculate penetration probability through the re- 
pulsive barrier, one of the terms we need for evaluation 
of fusion reaction rates , usin g the WKB approach and 
following [Bahcall et all (|1998[) . 

The fusion cross section of two isolated reacting nuclei is 
written as 



a{E) 



SjE) 
E 



-2 7rr)(B) 



(26) 



where S{E) is the astrophysical factor, E the center of 
mass energy of the fusing nuclei of charge Zi e and Z2 e 
colliding with relative velocity v/c = \f2EjiJu?), re- 
duced mass /i and 



^_ZlZ2e^ 

he 



1 



(27) 



/E V 2 

First of all we define the penetration factor when the pure 
Coulomb electrostatic potential energy barrier (r) = 
Zi Z2 (r) occurs, that is, when the reacting nuclei 
are isolated: 



: exp ' 



2 

he 



2^ic^ {v^{r)~E) 



1/2 



dr 



(28) 

where r^ is the classical turning point whose value is 
fixed by the relation (r^) = E. 

Secondly, by using the standard DH potential and still 
taking the turning point r-^-^ = r,, , a relation valid only 
for E > D, for small values of r^ /R^^ we have 



Tdh [E] = exp 



(29) 



If we consider the rates instead of the cross sections, fac- 
tor exp(— 7r?7r(-./i?j5jj) can be evaluated at most probable 
energy Eq in such a way that the rate can be factorized 
as the product of r^(£') times a factor fixed &i E — Ei^. 

Finally, we consider the deformed DH potential energy 
barrier 



Vq{r) 



D 



and the penetration factor 



1 - (1 - g) 



(30) 



exp ■ 



2 

he 



2/ic^ 



{Vq{r)-E) 



1/2 



dr 



(31) 

where must be derived from relation Vq{rq) — E. 

Because we consider q <1 potential energy Vq{r) has a 
cutoff and, as a consequence, < r < C,q/{1 — q). We 
write Eq. (|3ip as 



(32) 



r,(S) =cxp(-2^77(^)r, 



where function r^, which depends on quantity D / E, goes 
to one for D/E ^0. 

The evaluation of and Tq{E) can be worked out only 
numerically although, for small deformation (q ~ 1), the 
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Flc;. 1. — Log-linear jjlot of the quantity Tq/C^q as a function of 
D/E for several values of q (in the insert the region 10-1 < D/E < 
10 is expanded). 

penetration factor can be worked out analytically and 
is given as a product of (E) times a correction fac- 
tor. However we do not report here the expression for 
simplicity's sake. In the case q = 0, which represents 
the greatest deformation with respect to the exponential 
function, has the simplest analytical solution 



E 



Tq=0 



E + D 



(33) 



In Figure 1, we report the plot of quantity r^/Cg as a 
function of D/E for a few values of q between zero and 
one. 

In Figure 2, quantity 1 — is plotted as a function of 
D/E. 

Penetration of potential energy barrier Vq (r) at energy E 
is equivalent to penetration of the pure Coulomb barrier 
at an effective energy E + Uq where 



Uq = 



D 



rq{l/Kn) 



E 



(34) 



is a function of D/E. 

We have also calculated numerically equivalent energy 
Uq. The behavior of Uq/D^ as a function of D/E, for a 
few values of q is plotted in Figure 3. Finally, a quanti- 
tative comparison of quantities Tq/C,q, Uq/D and 1 — r^, 
corresponding to the values of q depicted in the figures 
can be obtained from the numerical Table 1. 

For q = {] potential Vq=Q is ctit at r = 2 (l/i^Dn)^^ 
we have Uq — D for any value of energy E. In the other 
cases of < g < 1 we have that: at high energy E > D, 
Uq approaches D. At low energy, E < D, Uq approaches 
the value (1 — q) D, meaning that Ui (DH exponential 
potential) approaches zero. Therefore, when the electron 
distribution is a deformed or generalized distribution of 
the g'-type, with a cut-off at r = (2 — g)/[(l — g) (l/R^^)], 
the penetration factor is enhanced even at very low en- 
ergy, except for the case g = 1. Enhancement depends 
on Uq which is an energy proportional to l/> 



Fig. 2. — Log-linear plot of the quantity 1 — as a function of 
D/E for several values of q (in the insert the region IQ-i < D/E < 
10 is expanded). 



1.0 



0.5- 



0.0 



-q = 0.00 
q = 0.25 
■q = 0.50 
■q = 0.75 
q = 0.99 



10"^ 



10" 



10^ 
D/E 



10=' 



Fig. 3. — Log-linear plot, in arbitrary unities, of Uq/D as a 
function oi D/E for several values of q. 



We have also interpolated the function Uq/D (for val- 
ues < g < 0.8) which is well described by the following 
analytic function 



g + g 



l-{l-q)a{D/Ey 



(35) 



where parameters a(g), h{q) and c(g) are given by 

a(g) = 0.484012 + 1.984961 g + 42.933001 g^ , 
6(g) = 0.958079 - 0.129717 g - 0.038993 g^ + 

0.804107 g'^ , (36) 
c(g) = 0.760938 - 1.272798 g + 0.493549 g^ . 
The interpolating function allows us to find the value of 
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TABLE 1 

Few numerical values of the quantities rq/(q, Uq/D and 1 — r,, in the energy range 0.1 < D/E < 30 for some values of 

<?■ 



D/E rq/Cq Uq/D 1 - Tq Tq/Cq Uq/D 1 ~ Tq Vq/Cq Uq/D 1 - Tq rq/(q Uq/D 1 - Tq 







q = 0.25 






q = 0.50 






q = 0.75 






q = 0.099 




0.10 


0.091 


0.988 


0.046 


0.091 


0.977 


0.045 


0.091 


0.966 


0.045 


0.091 


0.956 


0.045 


0.20 


0.167 


0.978 


0.085 


0.167 


0.958 


0.084 


0.168 


0.938 


0.083 


0.168 


0.920 


0.082 


0.30 


0.232 


0.969 


0.120 


0.233 


0.941 


0.118 


0.235 


0.915 


0.116 


0.236 


0.891 


0.114 


0.40 


0.288 


0.961 


0.151 


0.291 


0.927 


0.148 


0.294 


0.894 


0.145 


0.297 


0.866 


0.142 


0.50 


0.338 


0.955 


0.178 


0.343 


0.914 


0.174 


0.347 


0.877 


0.170 


0.351 


0.844 


0.166 


0.60 


0.382 


0.948 


0.203 


0.389 


0.902 


0.198 


0.395 


0.861 


0.193 


0.401 


0.825 


0.188 


0.70 


0.421 


0.942 


0.226 


0.430 


0.892 


0.220 


0.439 


0.846 


0.214 


0.447 


0.807 


0.208 


0.80 


0.457 


0.937 


0.246 


0.468 


0.882 


0.239 


0.479 


0.834 


0.232 


0.489 


0.792 


0.226 


0.90 


0.489 


0.932 


0.265 


0.503 


0.874 


0.257 


0.517 


0.822 


0.250 


0.529 


0.778 


0.243 


1.00 


0.518 


0.928 


0.283 


0.535 


0.866 


0.274 


0.552 


0.811 


0.266 


0.566 


0.765 


0.258 


2.00 


0.716 


0.896 


0.407 


0.763 


0.809 


0.393 


0.809 


0.735 


0.380 


0.850 


0.675 


0.367 


3.00 


0.826 


0.876 


0.482 


0.902 


0.774 


0.465 


0.977 


0.689 


0.449 


1.047 


0.621 


0.434 


4.00 


0.898 


0.863 


0.534 


1.000 


0.750 


0.516 


1.102 


0.657 


0.498 


1.198 


0.584 


0.482 


5.00 


0.949 


0.852 


0.572 


1.073 


0.731 


0.554 


1.200 


0.633 


0.535 


1.321 


0.556 


0.518 


6.00 


0.988 


0.844 


0.603 


1.131 


0.717 


0.583 


1.281 


0.613 


0.564 


1.426 


0.534 


0.546 


7.00 


1.019 


0.838 


0.627 


1.179 


0.705 


0.608 


1.349 


0.598 


0.588 


1.517 


0.516 


0.570 


8.00 


1.043 


0.832 


0.648 


1.219 


0.695 


0.628 


1.408 


0.584 


0.608 


1.597 


0.500 


0.590 


9.00 


1.064 


0.828 


0.665 


1.253 


0.686 


0.646 


1.461 


0.573 


0.626 


1.670 


0.487 


0.607 


10.00 


1.081 


0.824 


0.680 


1.283 


0.679 


0.661 


1.507 


0.563 


0.641 


1.735 


0.476 


0.622 


15.00 


1.140 


0.810 


0.733 


1.390 


0.652 


0.715 


1.684 


0.526 


0.695 


1.996 


0.434 


0.676 


20.00 


1.174 


0.801 


0.766 


1.459 


0.635 


0.749 


1.806 


0.503 


0.730 


2.188 


0.406 


0.711 



Uq{E) once we know the temperature and electron den- 
sity, fixed at Zi, Z2, and energy E for a certain value of 
q. In the range < D/E < 10^, for 0.8 < g < 1, func- 
tion (|35p still gives a good interpolation oiUq/D but with 
more complicated relationships of a(q), b{q) and c(q). We 
omit the details. 

Calculation of nuclear fusion rates (quantities weighted 
over the reacting nuclei distribution which in many cases 
is a generalized distribution with a proper ionic parame- 
ter qi close to one) requires the insertion into the average 
integral of the penetration factor that is a function of 
E. In the case of the pure Coulomb barrier the screen- 
ing factor can be factorized. In our general case this 
factorization is not possible and the behavior of Uq as 
a function oi D/E must be considered with care, both 
for non-resonant and resonant reactions. The same con- 
sideration is valid if instead of Uq we calculate the rates 
with the use of the plotted and tabulated function t{E). 



4. REPRESENTATIVE EXAMPLES OF PENETRATION 
PROBABILITY 

We report some representative examples on the eval- 
uation of electronic screening factor for non resonant 
and resonant fusion reactions of interest in solar core 
and in other dense astrophysical plasmas. Resonant fu- 
sions, beside electron scr eening, are influence d by a res- 
onance screening factor (jCussons et al.ll2002D . Further- 
more, Maxwellian rates can be corrected by non stan- 
dard ionic distributions. In this work, we are interested 
in modified DH potential and here we limit the discus- 
sion to electron screening factor. 

The enhancement of the penetration factor Tq{E) over 
the pure Coulomb penetration {E) can be expressed, 
using Eqs. ^ and dST]), by the ratio 



exp 



-2Tr 



^1^2 6^ /l 



^E+Uq ^/E 



(37) 

and the enhancement of VqiE) over T^^{E) = Tq^i{E) 

by 



fq, DH {E) - 



^q{E) 

^uAE) 



- exp 



\l^^{E)(Tq{E)-Tq^^{E))\ (38) 



It is evident from Figure 2 and from Table 1 that 
important enhancement over the DH potential comes 
from the energy range £^ ~ 2 • 10^"* 13 ^ i?, with a 
maximum at about E ~ 13/15. The most effective 
burning energy is at Eq = {Eq (fcr)^/4)^/^, where Eq 
is the Gamow energy. 

The first example concerns the p — p fusion 
in solar core. At the mass center energies above 
Eq ~ 5.9 • 10""^ MeV, Uq can be taken constant and 
equal to D = 0.7 • IQ-^MeV. The factor fq.dE) 
is of the order of few percent above unity at high 
energies, for ions with energy above Eq and belonging 
to the distribution tail and having fusion probability 
greater than those belonging to the head of distribution. 
However, at these high energies fq, {E) is practically 
equal to 1 also with q <^ 1. Therefore, the rate is not 
modified significantly respect to the standard value. At 
energies below Eq the enhancement is not negligible 
and Uq, which depends on E, goes to (1 — q) D as 
E ~* Q. We have, for instance, fq,^^{E — D) — 1.12 
and fq^^^{E = 0.1 D) = 1.65, when q = 0.5. However, 
protons with these energies have a very small probability 
to fuse and the rate in conclusion cannot change respect 
to the standard evaluation more than few percent, as 
also required by luminosity constraints. 
The above discussion is valid for all the other reactions 
of the hydrogen burning cycle, the effective burning 
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TABLE 2 

Penetration factors for ^^C - ^^C in laboratory experiments with graphite target. 



E(MeV) 










^q=0. 


,50 (il-') 


rq=0.99(-E) 


0.0059 


8.16 • 10" 


-494 


2.16 ■ 10" 


-349 


5.20 ■ 




1.75 ■ lO-'^'"'^' 


0.0590 


1.18 • 10- 


-156 


2.13 ■ 10" 


-149 


9.70 ■ 


10-141 


2.16 ■ 10-141 


0.5900 


4.91 • 10" 


-50 


8.63 • 10" 


-50 


2.00 ■ 


10-49 


1.97- 10-49 


2.2000 


2.92 ■ 10- 


-26 


3.16 • 10- 


-26 


3.56 ■ 


10-26 


3.16 ■ 10-26 



energy being in the range Eq ~ 15 30 KeV with 
D < Eo. 

Among the reactions of CNO cycle we consider ^'^N 
(p, 7) ^^O. The astrophysical factor rec ently measured 
(|Runkle et al.ll2005t llmbriani et al.ll2005f l has important 
consequences in the evolution of stars, estimation of age 
of globular clusters and of course in the evaluation of 
CNO neutrino flux (jLioliosI 120001 : iDegl'Innocenti et all 
[2OOI . At solar core conditions, we have that in the 
high energy region Uq = D <^ Eq ^ 27.0 KeV, therefore 
TqiE) « T^^{E). At low energy Uq ^ {1 - q) D and for 
E = D/15 = OMMeV we obtain Tq{E) = 2.63 • lO"* 
with fq,o^iE) = 3.33 at £; = 0.047 MeF and q = 0.5. 
Of course, /g, (-^) — > 1 in the g — + 1 limit as we can 
expect in laboratory experiments. 

Fusion reactions, like ^Li (p, a) a and ®Li (d, a) a are of 
particular interest for their implications in astrophysics. 
Experimental measurements in lab o ratory of the screen- 
ing potential (jEngstler et al.lll992t iPizzone et al.ll2003l ) 
have indicated a value of 350 400 eV, greater than 
the adiabatic theoretical one of 186 eV. We should 
evaluate D at the experimental conditions of target 
temperature and density. The quantity Uq is a fraction 
of D, decreasing to {1 — q) D and not to zero as in the 
standard DH approach. Assuming q = 0.5 we need 
D = 800 eV to obtain Uq ~ 400 e^. Unfortunately, 
we do not know the correct target temperature and 
density to calculate D. This difficulty is a motivation to 
consider questionable the application of this approach 
to laboratory experiments. 

We consider now the important non resonant reaction 
i^C (a, 7) 1^0 at a temperature of kT = 17.2 KeV 
and density p = lO'^'^ gr/cm'^ (as in the core plasma 
of helium burning red giant stars) with a mass 
fraction X(He) = 1/2 and A:(C) ^ 1/2. We have 
D = 0.573 KeV. With q = 0.5, Uq goes from about 
Uq = 0.286 KeV at very low energy to Uq = O.A96KeV 
at E ^ D and to Uq ~ D at high energy. The greatest 
value of fq. DH (E) is at E i ■ 10"^ KeV. At such 
low energy the screening factor /q.onl-^) ~ lO'^^, the 
Coulomb penetration factor being practically zero and 

r^E) « 10-458. 

Next we consider the example of "'^^C — ^^C fusion in 
laboratory experiments and in massive stars (in the 
classical thermonuclear regime, without considering, 
for simplicity's sake, effects du e to the presence of 
resonances and degeneratio n (jCussons et alJ l2002t 
lltoh et all 120031 : iFerrol l2004( l). In the experimental 
study of ^^C — ^^C fusion near the Gamow en- 
ergy, the target temperature is fc T ~ 6 • 10~^ MeV 
and the gra p hite density is p ~ 1.7 gr/ cm? 
(jSpillane et alJ l2007l ). therefore we have approxi- 
mately i?j3H ~ 0.31 • 10^ /to and D ~ 0.0168 AfeF. 



TABLE 3 

Penetration factors for ^^C - ^^C fusion in massive stars. 





i 9=0.25 (i-') 


i'9=0.99(-t'') 


0.0059 


1.89 ■ lO-"*^ 


3.07 ■ lO-'*^ 


0.0590 


7.50 ■ 10-41 


2.42 ■ 10-54 


0.5900 


2.19 ■ 10-^1 


1.08 ■ 10-^^ 


2.2000 


6.64 ■ 10-22 


2.62 ■ 10-22 



From (lAssenbaum et alJ [l987l) we have Uc = 5900 eV, 
i?a = 0.88 • 10V"i, rc (g) = exp(-87.21//E) and 
T^^{E) = exp(-87.21/\/E+T4). We assume q = 0.5 
representing an average deformation of electron dis- 
tribution. While Uc is fixed, Uq depends on energy 
E. Ratios T^„(E) / T ^(E) coincide with those of Ref. 
(jSpillane et al.ll2007f ). In Table 2 we report the values of 
r^(£;), T^^{E), Tq=o.Bo{E) and Tq=os9{E). Our value 
Tq^i{E) differs from V^^{E) because the first factor is 
calculated using Ug ^ijE) and the s econd using Uc fixed. 
In massive stars (jGasques et al.l 12005 ). the strong 
screening effect in a dense plasma can be simulated 
by assuming q = 0.25, a value that represents a large 
electron deformation. At fcT = 30 • 10~^MeF and 
p ~ 10'^ gr/ cm?, we have R^^ ~ 0.5 • 10^ /m and 
D ~ 1.04MeV. Enhancement of r,=o.25(-E) over 
rg=o.99(-E') is very evident; at -E = 2.2 MeV enhance- 
ment is given by a factor of about 2.5. In Table 3 we 
report, at several energies, Tq=o,25(E) and Fq^o. 99 (-£')■ 

Finally, in la supernova environmental conditions give 
kT ^ 5 - 10^^ MeV, central density p = 3 ■ 10^ gr/cm^ 
and X{C) = 1/3. Therefore, D = 1.5MeV, the greatest 
deviation is at -E = 0.1 MeV and consequently we have 
at this energy F^ = 7.9 • lO'^^s^ ^ 5 5 . ^0-78 ^^d 
Tq=o.5 = 4.5 • 10-^3 ^-^Yi fq,^^ = 6.9 • 10^. 
Oi course, in all above examples q is an arbitrary param- 
eter whose value should be determined a priori. Rates of 
fusion reactions will be evaluated using the parametric 
expression of Uq (E) given in Eq. ([55]) . 

5. CONCLUSIONS 

We have shown that in those systems where in the 
space coordinates stationary electron distributions devi- 
ate from the standard exponential one, the shielded elec- 
trostatic Coulomb potential is modified with respect to 
the standard DH potential derived using linear condi- 
tions and constraints. We used two different approaches 
that produced the same results. One consists of asso- 
ciating a Poisson equation having an unknown electron 
density distribution with a Bernoulli equation for quan- 
tity r Vq (r) whose solution is asymptotically a power law. 
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This result can also be obtained, and justified, by 
means of the super-statistics approach recently devel- 
oped within generalized statistical mechanics by consid- 
ering the inverse DH radius, as a fluctuating in- 
tensive parameter with relative variance to and parameter 
q given by relation q = 1-f w, that characterizes the fluc- 
tuation. 

The g-modifled DH potential has asymptotic power law 
behavior and we have discussed its meaning for the case 
q < \ which holds when the electrostatic potential has 
a depleted and cut tail. Fluctuation of the inverse DH 
radius may be due to temperature and density fluctua- 
tion of the electrons surrounding the reacting nuclei in 
astrophysical plasmas. 

Exact evaluation of the penetration factor has to be car- 
ried out numerically. However, in the case of small defor- 
mations, one could calculate T{E) as the pure Coulomb 
penetration factor (E) times a correction that may 
sensibly differ from the standard DH correction. 
We have reported detailed plots and numerical tables 
of several useful quantities that allow evaluation of the 
penetration factor and the reaction rates which, being 
weighted integrals, depend sensibly on the behavior of 
the penetration factor as a function of E. 
Penetration factor rq{E) can be given as the pure 
Coulomb penetration factor at an equivalent energy 
E + Uq. Energy Uq is not a constant and can be eval- 
uated only numerically. We have given a useful fit of 
it. Uq is a function of energy D, characteristic of re- 



acting nuclei, their density and fusion temperature, the 
constant of proportionality depending on q and varying 
with D/E. Uq is proportional to l/V/cT and to ^/n as 
observed experimentally in fusion reactions in metal ma- 
trices. The complete and correct expression of the pene- 
tration factor or of Uq or 1 — Tq reported in the paper are 
necessary to evaluate the fusion rates at any energy. In 
fact, a possible deformed ion distribution together with 
an electron deformed distribution may sensibly affect the 
values of the rates. In the evaluations of the rates of fu- 
sion reactions, such as in the examples of p — p, ^^N 
(p, 7) 1^0, ^Li (p, a) a and ^Li (d, a) a, (a, 7) 
^^O and ^^C — ^^C fusion, mentioned in the previous 
paragraph, important enhancements may occur, in addi- 
tion to possible enhancements or decrements due to non- 
Maxwellian energy-momentum ion distributions, if the 
electron clouds surrounding the reacting nuclei are spa- 
tially modified by an assumed q < 1 distribution. Spe- 
cific applications of our approach are in progress; they are 
comprehensive of discussion of the recent experimental 
results of the LUNA collaboration (,Spillane et al.. .200'^ : 
llmbriani et al.l [20051 : lGvurkvll2007D to derive astrophys- 
ical factors at stellar energy range and will be the ar- 
gument of a further work. We are confident that our 
approach and numerical tables shown here will be used 
to study astrophysical process where temperature and 
density fluctuations cannot be neglected to evaluate DH 
shielding effect on reaction rates. 
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